arXiv:astro-ph/9706053vl 5 Jun 1997 


Astrophysical Sources of Stochastic Gravitational Radiation 

in the Universe 


K.A.Postnov 

Sternberg Astronomical Institute, Moscow State University, 

Moscow 119899, Russia 

Abstract 

Stochastic gravitational waves (GW) associated with unresolved astrophysi¬ 
cal sources at frequency bands of the ongoing GW interferometers LIGO/VIRGO 
and LISA are studied. We show that GW noise from rotating galactic neutron 
stars with low magnetic fields may reach the advanced LIGO sensitivity level 
at frequency / ~ 100 Hz. Within LISA frequency band (10“^ — 10“^ Hz), the 
GW background from galactic binary stars is shown to mainly contribute up 
to a frequency of 3 x 10“^ Hz, depending on the galactic rate of binary white 
dwarf mergers. To be detectable by LISA, relic GW backgrounds should be as 
high as Hgw^ioo > 10“^ at 10“^ Hz. 


1 Introduction 

In a few years, with the completion of construction of gravitational wave detectors 
of high sensitivity, a new window into the Universe will be open (see Thorne 1995, 
Schutz 1996 for a recent review of gravitaitonal wave astronomy). In this connec¬ 
tion, extensive studies of possible sources of gravitationa radiation are now being 
conducted. The most promising targets for the initial laser interferometers with the 
rms sensitivity level hrms ~ 10“^^ at / = 100 Hz are coalescing binary neutron stars 
and/or black holes, which can be observable from distances up to 100 Mpc (Lipunov, 
Postnov & Prokhorov 1997). 

Stars are most numerous baryonic objects in the Universe (~ 10^^ within the 
Galaxy, ~ 10^^ within the Hubble radius Rh SOOOh^oo Mpc, where hioo = Hq/100 
km/s/Mpc is the present value of the Hubble constant). Among them the most 
signihcant sources of GW are rotating triaxial neutron stars (NS) and binary stars. 
Only a small fraction of galactic NS (about 700) is observed as radiopulsars, so when 
trying to search for GW from them we have an advantage of knowing the precise spin 
period and the position on the sky. The same argument relates to the binary stars 
with known orbital periods. But the vast majority of them are unresolved sources 
and will form a stochastic background. A stochastic GW background is commonly 
measured in terms of the energy density per logarithmic frequency interval related 
to the critical energy density to close the Universe, Utcw = dEcw/dhr f /PcrC^ {pcr = 
SHq/StiG ~ 1.9 X 10“^®h^QQ g cm“^ where c is the speed of light). For comparison with 
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dimensionless detector’s sensitivity h, one commonly uses the equivalent characteristic 
strain h^if) = (l/27r)(ifo//)^dJ/^. 

Being interesting by themselves, astrophysical backgrounds, however, are viewed 
as a noise burying a possible cosmological gravitational wave background (CGWB), 
which bears the unique imprint of physical processes occurring at the very early 
(near-Plankian) age of the Universe (see e.g. Grishchuk 1988 for a review). Exact 
value of GGWB is still very controversial (see Grishchuk 1996 for fresh estimates). 
What is more reliable (however, not completely parameter-free) are GWBs formed by 
known astrophysical sources (old neutron stars, binary stars), and here we address the 
question how much astrophysical GW noises contribute at LIGO/VIRGO (10-1000 
Hz) and LISA (10“"^ — 10“^ Hz) frequencies. 

2 GW noise from sources with changing frequency 

To calculate GW noise produced by some unresolved sources we need to know the 
number of sources per logarithmic frequency interval. At the hrst glance, this would 
require knowing precise formation and evolution of sources. However, when only GW 
carries away angular momentum from the emitting source, the problem becomes very 
simple and physically clear. 

GW energy loss leads to changing the frequency u = 2nf of emitting objects. 
In the case of a rotating triaxial body, the positive rotational energy Erot = 

(here I is the moment of inertia) is being lost and the spin frequency (hence, GW 
freqency) decreases. In contrast, in the case of a binary star with masses of the 
components Mi AND m 2 in a circular orbit ofr adius a the negative orbital energy, 
Eorh = -MiM2/2a ~ (here M = {Mi + M 2 f/\MiM 2 f/^ is the so-called 

“chirp mass” of the system) is being lost and the orbital frequency increases. 

To a very good approximation, the conditions of star formation and evolution of 
astrophysical objects in our Galaxy may be viewed as stationary. This is true at least 
for last 5 billion years. Let the formation rate of GW sources be TZ. For example, 
the mean formation rate of massive stars (>10 Mq to produce NS) is about 1 per 30 
years. 

The stationarity implies that the number of sources per unit logarithmic frequency 
interval is 

dN/d\nf^N{f)=nx{f/f). (1) 

The total energy emitted in GW per second per unit logarithmic frequency interval 
at / by all such sources in the galaxy is 

dEcw / {dt din f) = L(/)gw = L{f)GwN{f) = L{f)Gw'^ x (///) , (2) 

where L{f)GW is the GW luminosity of the typical source at frequency / {L{f)Gw oc 

for non-axisymmetric neutron stars, L(/)gw oc for binary systems). 
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Finally, for an isotropic background we have 

^Gw{f)PcrC^ = ^(/)Gw/(47rc(r)2) (3) 

where (r) is the inverse-square average distance to the typical source. Strictly speak¬ 
ing, this distance (as well as the binary chirp mass Ai and moment of inertia I of 
NS) may be a function of frequency since the binaries characterized by different At 
may be differently distributed in the galaxy. We are highly ignorant about the real 
distribution of old NS and binaries in the galaxy, but taking the mean photometric 
distance for a spheroidal distribution in the form dN oc exp[—r/ro] eyi^[—{z/Z qY] (r 
is the radial distance to the galactic center and 2 : the hight above the galactic plane) 
with ro = 5 kpc and Zq = 4.2 kpc with (r) ^ 7.89 kpc is sufficient for our purposes. 

For the cases considered the energy reservoir radiated in GW is either rotational 
energy (neutron stars) or orbital energy (binary systems), both depending as some 
power of the corresponding frequency: E oc /“, ajvs = 2, a/,s = 2/3. Hence the 
frequecy change (///) may be found from the equation dE/dt = aE{f/f) By energy 
conservation law dE/dt = {dE/dt)cw + {dE/dt)EM + {dE/dt),,, where index EM 
stands for elecromagnetic losses and ... means other possible losses of energy. Finally, 
we obtain 

(///) = aE/{{dE/dt)Gw + {dE/dt)EM + {dE/dt)„) (4) 

and ^ 

L{f)GW = —(dE/dt) em I (dE/dt),,, 

^ {dE/dt)GW (dE/dt)GW 

The remarkable result is that if GW is the dominant source of energy removal, the 
resulting GW stochastic background depends only on the source formation rate: 

^Gw{f)PcrC^ = TZaE / {incir)"^) (6) 

3 GWB from old neutron stars at LIGO/VIRGO 
frequencies 

Spin evolution of rotating non-axisymmetric NS with ellipticity e and magnetic mo¬ 
ment fi may be driven by GW or elecgtromagnetic losses. The condition that a 
stochastic signal appears within the detector band deoends on the rate of the fre¬ 
quency change. The upper frequency of the stochastic background for pure electro¬ 
magnetic energy losses is ~ 10^(Hz)773o^/|5^/i7o^ where /iso = /i/(10^°G cm^) is 
NS magnetic moment. For pure GW losses this upper frequency is ^ 1.4 x 
10"‘(Hz)773o^/45^'^"^e_e'^^ where e_6 = e/10® is the NS ellipticity. 

For plausible values of the NS magnetic helds (/iso = 10“^-10^) and ellipticities 
(e_6 = 10“®-10^), at any frequency < 10® Hz we deal with stochastic backgrounds 
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from galactic NS. Physically, this is due to the inability of old NS to leave frequency 
interval Ao; ~ lo during the typical time between consecutive supernova explosions. 

For purely GR-driven NS spin-down the resulting spectrum is independent of the 
unknown value of e in the NS population. Any additional braking mechanism always 
lowers the resulting signal. Taking typical values J = 10^^ g cm^, 'R. = 1/30 yr“^ we 
obtain from Eq. 

Gjvs ~ 10-^7^^/'/45(//100i7^)2h^oo(r/10A;pc)-2 (7) 

or in terms of he 

K « jsfaT^ « 10 -"- ( 8 ) 

(here we assumed the characteristic distance to NS population of order 10 kpc). 
Remarkably, this limit does not depend on frequency. The GR background of such 
strength could be detected by the advanced LIGO/VIRGO interferometers in one 
year integration (Thorne 1987; Giazotto 1997). 

For realistic NS parameters the ratio of electromagnetic to GW losses x = Eem/Egw 
is 

P o /l00Hz\^ 

X ~ 4000/i3Qe_g f j j (9) 

and electormagnetic losses becomes insigniheant (x -C 1 ) only at high frequencies 
f > fer ^ 6.3(kHz) ^ If we would take e _6 = 10“^ and /X 30 = 10“^ as in millisecond 
pulsars, we would obtain f^r ~ 630 Hz, however millisecond pulsars are spun up by 
accretion in binary systems and are not considered here. 

Therefore, for realistic NS we must consider the case x 3> 1. Then the stochastic 
background from old NS becomes 

hAj) ~ 5 X 10-=* (10) 

and is lis below even advanced LIGO sensitivity at / ~ 100 Hz. 

We have shown that if the NS form ellipticity is present, the stochastic GR back¬ 
ground produced by old NS population is naturally formed due to NS rotation brak¬ 
ing. In the limiting case when only GR angular momentum loss causes NS spin-down, 
this background is independent on both exact value of the NS form ellipticity e and 
frequency and can be detected by advanced LIGO/VIRGO interferometers. In re¬ 
ality, the magnetic held of NS causes more effective electromagnetic energy loss: 
to be insigniheant, the magnetic held of a NS should be less than (see Eq. (i) 

/i < 1.5 X 10^®(G cm^)e_6/ 

According to Urpin & Muslimov (1992), the magnetic held can decay very fastly 
provided that the held was initially concentrated in the outer crust layers with the 
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density < 10^° — 10^^ g cm“^, and such very low magnetic field for old NS may be 
possible. In the limiting case that the NS magnetic field does not decay at all (for ex¬ 
ample, if only accretion-induced field decay is possible in binary systems (Bisnovatyi- 
Kogan & Komberg 1974)), old NS should lose their energy through electromagnetic 
losses and be very slow rotators with periods of about a few seconds. Then the ini¬ 
tial magnetic field distribution becomes crucial. If it is centered at 10^2 G (as 
implied by radipulsar P P measurements), we have little chances to detect the old 
NS population at 10-100 Hz frequency band unless close mean distances (<10 kpc) 
are assumed (Giazotto et ah 1997). However, if nature prefers a scale-free law (like 
f{fi) oc l//i), the fraction of low-field NS could amount to a few 10% and they can 
contribute to the frequency-independent GR background. Then Eq. (|^) implies that 
such a background can be detected by the advanced LIGO/VIRGO interferometer 
in the frequency band 10-1000 Hz in one-year integration even if the formation rate 
of such NS is as small as 1 per 300 years and the characteristic distance to them is 
100 kpc. 


4 GW noise from unresolved binary stars at LISA 
frequencies 

Inside LISA frequency range, 10“^ — 10“^ Hz, only coalescing binary white dwarfs 
(WD) and binary neutron stars contribute. Even if binary neutron stars coalesce at 
a rate of 1/10000 yr in the Galaxy, their number still should be much smaller than 
the white dwarf binaries, and in this section we restrict ourselves to considering only 
binary WD. 

Substituting E = Eorb ~ into equation (H) we obtain 

^wd{!) ~ 2 X 10-S77ioo(//10-'Hz)2/3(Wl/MQ)5/3^(r)/10kpc)-2h7o2 , (n) 

where 77ioo = 77/(0.01 yr“^) is the galactic rate of binary WD mergers. 

In terms of the characteristic dimensionless amplitude of the noise background that 
determines the signal-to-noise ratio when cross-correlating outputs of two independent 
interferometers we have 

K{f) ^ 7.5 X 10-2°77;g(//10-3Hz)-2/3(A:i/Mo)''/®((r)/10kpc)-i (12) 

Equation (|l^ shows that at high frequencies of interest here the GW background 
is fully determined by the galactic rate of binary WD mergers and is independent 
of (complicated) details of binary evolution at lower frequencies (the examples of 
calculated spectra at lower frequencies see in Lipunov & Postnov 1987; Lipunov, 
Postnov & Prokhorov 1987; Hils et ah 1990). 

But the real galactic merger rate of close binary WD is unknown. One possible 
way to recover it is searching for close white dwarf binaries. A recent study (Marsh 
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et al. 1995), revealed a larger fraction of such systems than had previously been 
thought. Still, the statistics of such binaries in the Galaxy remains very poor. 

If coalescing binary WD are associated with SN la explosions, as proposed by Iben 
& Tutukov (1984) and further investigated by many authors (for a recent review of 
SN la progenitors see Branch et al. 1995), their coalescence rate can be constrained 
using much more representative SN la statistics. Branch et al. (1995) concluded 
that coalescing CO-CO binary WD remain the most plausible candidates mostly 
contributing to the SN la explosions. The galactic rate of SN la is estimated 4 x 10“^ 
per year (Tamman et al. 1994; van den Bergh and McClure 1994), which is close to 
the calculated rate of CO-CO coalescences (~ (1 — 3) x 10“^). The coalescence rate 
for He-CO WD and He-He WD (other possible progenitors of SN la) falls ten times 
short of that for CO-CO WD (Branch et al. 1995). As SN la explosions may well be 
triggered by other mechanisms, we conclude that the observed SN la rate provides 
a secure upper limit to the double WD merger rate regardless of the evolutionary 
considerations. 

The upper limit ( [T^) is plotted in Fig. 1 for different rates of binary WD mergers 
T^-ioo = 1,1/3,1/10,1/30 assuming the chirp mass M. ~ O.52M0 (as for two CO white 
dwarfs with equal masses Mi = M 2 = O.OMq). These lines intersect the proposed 
LISA rms sensitivity at / > fum ~ 0.03 — 0.07Hz. This means that at frequencies 
higher than 0.07 Hz no continuous GW backgrounds of galactic origin are presently 
known to contribute above the rms-level of LISA space laser interferometer. The 
contribution from extragalactic binaries is still lower regardless of the poorly known 
binary WD merging rate (at least in the limit of no strong source evolution with 
z). Other possible sources could be extragalactic massive BH binary systems (e.g. 
Hils and Bender 1995). Their number in the Universe can be fairly high (e.g. Rees 
1997), but no reliable estimates of their contribution are available at present. The 
lower limit (§) is already close to the LISA sensitivity limit at 0.1 Hz, but we stress 
that the assumptions used in its derivation are upper limits, so the actual frequency 
beyond which no binary stochastic backgrounds contribute may be three times lower. 
This precise limit depends on the details of binary WD formation and evolution which 
are still poorly known. 

Fig. 1 demonstrates that the calculated GW background intersects LISA sensitiv¬ 
ity curve at frequencies ~ 0.05 Hz, and Bender et al’s curve at even lower frequencies 
~ 0.01 Hz. The latter is probably due to Bender et al’s curve being derived from 
observational estimate of double WD galactic density in the solar neighborhoods; we 
stress once more that once formed, the binary WD will evolve until the less massive 
companion hlls its Roche lobe; unless the mass ratio is sufficiently far from one (cf. 
Webbink 1984), the merger should occur. Therefore Bender et al’s curve provides a 
secure lower limit to the galactic binary stochastic GW background. 

Presently, we cannot rule out the high galactic double WD merger rate (1/300- 
1/1000 yr“^), and therefore can consider fum to lie within the frequency range 
0.01 — 0.07 Hz. We conclude that no GW background of galactic origin above this 
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Figure 1: Galactic binary GW background he as given by Bender (1996) (filled 
quadrangles) and calculated for a model spiral galaxy with the total stellar mass 10 ^^ 
Mq (the solid curve). Average photometric distance 7.9 kpc is assumed. Thin straight 
lines marked with 100, 300, 1000, 3000 are the analytical upper limit (eq. [0) for 
binary WD merger rates 1/100, 1/300, 1/1000, and 1/3000 yr“^ in a model spiral 
galaxy, respectively, assuming M. = 0.52 Mq. Straight dashed lines labeled by 10“®, 
10“® show GW backgrounds corresponding to constant Vtcw- The proposed LISA 
rms noise level (hrms) and sensitivity to bursts hsB = hy/bhrms are also reproduced 
(cf. Thorne 1995; Fig. 14). 


frequencies should contribute at the rms-noise level of LISA interferometer, and hence 
the detection of an isotropic stochastic signal at frequencies 0.03 — 0.1 Hz with an ap¬ 
preciable signal-to-noise level (which possibly may be done using one interferometer) 
would strongly indicate its cosmological origin. To be detectable by LISA, the power 
of relic GW background should be Ggiv^ioo > this frequency range. 
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